Comment on: Role of Intermittency in Urban 
Development: A Model of Large-Scale City For- 
mation 

In Ref. |jj a model for large-scale city formation was pro- 
posed. It is based on the discreet dynamics (eqs. (1,2) of 
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n i (t + l) = (l-a)n i (f)/ i (t) + ! £ n 3 (t)f 3 (t), (1) 
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where rn(t) is the population at site i of a square lattice, 
and fi(t) > are random multiplication factors drawn 
independently from a distribution P(f) with the property 
(/) = 1. The authors claim that such a multiplicative 
process with diffusion gives rise to a stationary state with 
a power law distribution of city sizes n, P(n) ~ n~ r , 
where r ~ 2 in their simulations is in agreement with 
empirical observations or real city size distribution. 

We demonstrate here that i) Eq. (|l|), with (/) = 1, 
does not lead to a stationary state; ii) the probability 
distribution of rij does not have a power-law tail; Hi) both 
stationarity and a power-law tails of the distribution are 
recovered if a constant source term is added to the RHS 
of the Eq. (Q), but then the uniqueness of the exponent 
t = 2 is lost. 

We have performed a numerical simulation of the dy- 
namical rules (^) and found that, contrary to the claims 
of Ref. [H , these dynamical rules do not lead to a station- 
ary state. Instead, the total population N(t) = 
decays to zero exponentially in time. For example, for 
a = 0.25, and / being equal to 2 or with equal prob- 
abilities, we found that after 1000 time steps N(t) is of 
order 10~ 15 in clear contradiction with the results of Ref. 
|jj for the same set of parameters. This is not a finite size 
effect (we found the same rate of decay for sizes in the 
range 64 2 -200 2 ). This inconsistency suggests that there 
is some ingredient in the model simulated by Zanette and 
Manrubia, which was not reported in their paper. 

To understand better the reasons for this exponential 
decay let us consider first the case a = 0, where the 
population rij(t) at each site undergoes a random mul- 
tiplicative process Uj(t + 1) = ni{t)f%(t) and is uncou- 
pled from other sites. For (/) = 1 the expectation value 
(rii(t)) does not change in time. It is known, however, 
that in such process, (rii(t)) is dominated by extremely 
rare events, when rn is exponentially large in t. By the 
virtue of Central Limit Theorem, in any typical realiza- 
tion, ni(t) ~ e^ 11 -^*. One can show that for any P(f) 
such that (/) = 1, one has (In/) < 0, so that the typical 
rii(t) vanishes exponentially for t — > oo. The diffusion 
(a > 0) alone cannot reverse this typical decay. Indeed, 
the total population also undergoes a multiplicative pro- 
cess: N(t + 1) = F(t)N(t), where F = J2ifi n i/ N is 
the "population" average of /j. Again (F) = 1 implies 
that (InF) < 0, so that typically N(t) ~ e {lnF)t -> for 
t — » oo. 

Note that Eq. ([!]) is the equation for the partition 
function J2i n i(t) °f a directed polymer of length t in a 


three-dimensional random media Q, and — (InF) is the 
polymer's free energy per unit length. For this problem 
it is known [|| that P(n) is not a power law but rather a 
log-stretched-exponential law. 

On the other hand, if a positive constant is added to 
the right hand side of Eq. (|l|) , it describes a multiplica- 
tive noise process with a lower wall. The purpose of this 
extra source term is to prevent ni(t) from becoming too 
small, while for large its influence can be safely ne- 
glected. Such rules dynamics have received much at- 
tention recently both for a > ||, and a = [Q. 
For a = it was shown Q that a stationary state ex- 
ists, provided (In /) < 0. The distribution of n has a 
power law tail n~~ T with r determined by (/ T_1 ) = 1. 
Clearly (/) = 1 implies r = 2, but the condition (/) = 1 
is no longer necessary for the existence of a stationary 
state, and cannot be justified on these grounds. For 
a > wc found that the system is stationary (it is in 
the "pinned" phase ||) and P(n) still has a power law 
tail with t = 2. This can be justified by the fact that the 
correlation length £(a) is finite away from the depinning 
transition ||. Therefore, a large system can be divided 
into many virtually uncoupled blocks, whose total pop- 
ulation fi x — J2\i- X \<£(a) n i undergoes a multiplicative 
dynamics n x (t + 1) = f x (t)n x (t) with (/} = 1. The last 
condition ensures that r = 2 for any a. If (/) ^ 1 instead 
r depends on a. 

In short, the results of ref. are inconsistent with 
their rules of dynamics (Q), but compatible with the be- 
havior of models with a lower wall Q . For such models a 
stationary state does not require (/} = 1, and any power 
law exponent r > 1 is feasible. 
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